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Abstract — In  this  paper,  we  address  the  problem  of  find¬ 
ing  concertina  wire  in  three-dimensional  (3-D)  data.  Wire 
entanglements  constitute  a  major  obstacle  to  the  mobility  of 
Unmanned  Ground  Vehicle  because  of  their  widespread  use 
and  the  difficulty  to  detect  them.  We  pose  the  problem  in  term 
of  finding  thin  structures  organized  in  complex  patterns.  Such 
problem  did  not  received  as  much  attention  as  linear  and 
planar  structures  segmentation.  We  are  interested  especially 
in  the  problems  posed  by  repetitive  and  symmetric  structures 
acquired  with  a  laser  range  finder.  The  method  relies  on  3-D 
data  projections  along  specific  directions  and  2-D  histograms 
comparison.  The  sensitivity  of  the  classification  algorithm 
to  the  parameter  settings  is  evaluated  and  a  segmentation 
method  proposed. 

I.  Introduction 

Wire  entanglements  are  commonly  used  to  reinforce 
natural  obstacles  and  to  impede  forces  in  military  oper¬ 
ations  both  in  urban  and  natural  environment.  Concertina 
wires  are  such  an  obstacle.  Typically,  it  is  made  of  a  thin 
cable  (3.5  mm  diameter)  forming  a  one  meter  diameter 
coil.  It  can  be  fitted  with  blades  (10  mm  x  13  mm)  at 
regular  interval  along  the  cable  (every  2.6  cm).  Coils  can 
be  arranged  in  various  ways:  supported  by  poles,  stacked 
on  top  of  each  other  to  create  a  fence,  or  just  laid  across 
a  street.  As  such,  concertina  wires  constitute  very  difficult 
but  omnipresent  obstacles  to  detect  by  Unmanned  Ground 
Vehicles  (UGV). 

In  this  paper,  we  investigate  the  detection  of  concertina 
wires  in  three-dimensional  (3-D)  laser  data.  Effective  meth¬ 
ods  have  been  developed  for  classification  and  segmenta¬ 
tion  of  natural  scenes  from  such  data.  Applications  are 
scene  interpretation,  obstacle  detection  or  terrain-based 
localization  [1],  [2],  [3],  [4].  A  substantial  body  of  work 
exists  also  on  how  to  extract  planar  and  linear  objects 
from  scattered  3-D  point  clouds,  see  for  example  [5], 
[6].  Methods  were  even  proposed  to  detect  solid  surfaces 
hidden  inside  foliage  [7].  But  small  scale  structures,  such  as 
wires  and  thin  poles,  are  still  challenging  for  segmentation 
algorithms. 

The  problem  of  concertina  wire  detection  is  an  instance 
of  a  larger  class  of  problems  involving  the  segmentation 
of  and  analysis  for  shape  retrieval  in  3-D  points  clouds. 
This  problem  is  particularly  acute  with  porous  scenes 
(vegetation)  that  include  in  addition  thin  structural  elements 
(branches,  poles).  In  this  paper,  we  look  at  the  detection 
and  segmentation  of  organized  thin  structures.  In  our 
terminology,  thin  structures  differ  from  linear  structures  in 


the  following  ways:  linear  structures  are  elongated,  locally 
smooth  set  of  points  (typically  large  branches  or  posts) 
while  thin  structures  are  elongated  and  a  few  pixels  wide 
(typically  small  branches  or  wires).  Organized  thin  struc¬ 
tures  are  complex  patterns  of  thin  structures.  There  are  two 
key  issues  in  detecting  such  structures:  1)  the  small  number 
of  points  on  the  target  of  interest,  and  2)  their  structural 
similarity  with  vegetation  and  foliage.  We  are  interested  in 
characterizing  such  structure  by  building  a  discriminative 
local  structure  signature.  We  focus  our  effort  on  a  special 
class  of  organized  thin  structures:  periodic  and  symmetric 
structures. 

Figure  1  shows  a  representative  scene  containing  foliage, 
grass,  a  post  and  a  concertina  wire. 


(a)  (b) 


Fig.  1.  (a)  Representative  terrain  of  interest  containing  vegetation  and 

a  concertina  wire,  (b)  Close-up  view  of  a  concertina  wire. 

Our  proposed  local  structure  signature  relies  on  data 
reduction,  both  in  dimension  and  in  size,  by  projecting 
and  accumulating  3-D  data  into  different  2-D  sub- spaces. 
The  approach  is  inspired  by  the  work  on  local  shape 
representation  with  Spin-Images  by  Johnson  [8]  and  with 
Shape  Contexts  by  Belongie  [9].  The  level  of  organization 
of  the  point  cloud  in  which  we  are  interested  is  quantified 
by  an  absolute  measure  of  local  bilateral  symmetry. 

The  rest  of  the  paper  is  composed  of  four  sections.  In 
the  next  section  we  review  current  methods  on  sensing  and 
segmentation  of  thin  structures.  In  Section  III,  we  present 
the  details  of  our  approach.  In  Section  IV,  we  illustrate  the 
performances  of  the  method  on  real  data  before  concluding. 

II.  Related  work 

A.  Sensing 

The  two  closest  related  fields  in  thin  structure  detection 
is  power  line  detection  for  aircraft  obstacle  avoidance  at 
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Fig.  2.  FITG  field  demonstration,  (a)  Robot  in  a  forest  area  moving  toward  a  one  inch  wire  hang  between  two  trees,  (b)  A  close-up  view  of  the  wire, 
(c)  Classification  results.  In  red/blue/green  3-D  points  classified  as  surface/linear/vegetation. 


low  altitude  and  stand-off  trip  wire  detection  for  land¬ 
mine  clearance.  Laser-based  obstacle  detection  systems  are 
commercially  available  for  power  line  detection  on-board 
manned  helicopters.  Typically,  they  are  designed  to  detect 
at  600  m  cables  with  diameter  of  several  centimeters, 
stretching  in  a  linear  fashion  between  towers  with  no 
clutter.  The  emphasis  is  on  early  and  reliable  detection  and 
not  on  3-D  reconstruction  [10].  Other  sensing  modalities 
tested  include  millimeter  wave  radar  [11]  and  passive 
infrared  camera  [12].  Trip-wires,  used  in  some  landmines, 
can  be  detected,  using  passive  polarimetric  infrared  sensors 

[13]  or  active  illumination  with  short-wave  infrared  laser 

[14] .  Those  wires  are  a  few  millimeters  in  diameter  and 
cluttered  by  vegetation. 

The  accurate  3-D  reconstruction  of  thin  structures  is 
very  difficult  in  natural  cluttered  settings  such  as  the  ones 
encountered  by  ground  mobile  robots.  Thin  structures  will 
produce  mixed  pixels  that  are  difficult  to  compensate  for 
using  methods  based  on  local  smoothness  of  the  scene  [15]. 
Tuley  in  [16]  reported  interesting  results  on  analysis  and 
removal  of  artifacts  in  3-D  LADAR  data  for  simple  natural 
scenes. 

Only  one  publication  dedicated  to  the  detection  of  con¬ 
certina  wires  could  be  found.  In  [17],  the  authors  present  a 
device  to  be  effective  at  imaging  glow-contrast  underwater 
targets. 

B.  Segmentation 

In  this  paper,  we  address  the  problem  of  structure 
detection  and  segmentation  from  3-D  point  clouds  collected 
from  a  single  sensor  location  or  integrated  from  multiple 
locations.  In  [2],  Lacaze  proposed  to  use  the  permeability 
of  the  scene  to  laser  beam  to  discriminate  between  solid 
surfaces  (ground,  obstacle)  and  vegetation.  This  method  is 
popular  for  outdoor  robot  navigation  [18],  [1],  [19]  but 
does  not  allow  object  segmentation  or  scene  structures 
inference.  Methods  that  fit  primitives  to  point  clouds  are 
difficult  to  use  practically  for  large  data  sets  containing 
multiple  complex  structures,  in  opposition  to  multiple 
planar  surfaces  that  can  be  handled.  Similar  methods  using 
a  single  range  image  are  sensitive  to  occlusion.  In  addition 
they  focus  on  piece- wise  smooth  structures,  not  on  porous 
objects  such  as  a  concertina  wire.  Tensor  voting  was  proven 
to  be  a  powerful  tool  to  segment  surface,  curve  and  junction 
from  3-D  data  [6].  We  follow  a  similar  approach  in  [20], 


as  described  in  the  next  section.  Tubular  structures  from 
veins  can  be  extracted  from  dense  3-D  volumetric  medical 
data  [21].  Data  that  differs  very  much  from  the  sparse 
discontinuous  data  we  work  with. 


C.  Our  previous  approach 

In  [20],  we  proposed  a  method  to  classify  3-D  LADAR 
data  in  natural  terrain  into  three  classes:  vegetation,  solid 
surface  and  linear  structure.  The  method  estimates  the  local 
point  distribution  in  space  and  uses  a  Bayes  classifier  to 
produce  the  probability  of  belonging  to  each  class.  Priors 
are  modeled  as  Mixtures  of  Gaussians  and  parameters 
are  learned  using  the  Expectation-Maximization  algorithm 
(EM)  on  labelled  data.  At  each  point,  the  scatter  matrix  is 
computed  using  a  predefined  support  region.  The  principal 
components  of  this  matrix  are  used  to  define  three  saliency 
features  for  each  scale  [6],  characterizing  the  3-D  points’ 
spatial  distribution  into  the  three  classes. 

We  demonstrated  this  approach  for  general  purpose 
terrain  classification  and  wire  detection  with  the  GDRS 
XUV  unmanned  vehicle.  Tests  were  conducted  at  Fort 
Indiantown  Gap  in  the  spring  of  2004.  Figure  2  show 
a  typical  result  obtained.  The  processing  was  done  on¬ 
board  of  the  XUV.  The  XUV  was  teleoperated  and  directed 
toward  a  wire  at  roughly  2  m/s.  Results  of  the  on-board 
processing  was  send  back  to  the  operator  control  unit  via 
a  wireless  link.  The  data  exchange  between  the  different 
processes  (laser  data  acquisition,  processing,  display)  was 
done  though  shared  memory  using  the  Neutral  Message 
Language  (NML). 

This  method  gives  good  results  with  thin  linear  objects 
with  low  curvature.  Multiple  thin  linear  objects  close  to 
each  other,  intersecting  or  forming  junctions  are  more 
challenging.  Thin  structures  forming  a  compact  pattern 
such  as  a  coil  fall  into  that  category.  Figure  9-(c)  shows  an 
example  where  a  complex  symmetric  structure  cannot  be 
detected  reliably  with  such  a  method.  The  scene  contains  a 
concertina  wire  cluttered  by  vegetation.  Most  of  the  points 
from  the  structure  are  classified  as  scatter  (green  points) 
while  few  are  as  linear  structure  (blue  points)  as  we  would 
like  to  see.  In  the  next  section  we  present  a  method  to 
address  this  problem. 
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Bilateral  symmetry 


Fig.  3.  (a)  Projection  for  bilateral  symmetry  detection,  (b)  Example  of  data  from  a  concertina  wire  projected  into  two  plane  configurations:  the 

top/bottom  2-D  histograms  correspond  to  the  plane  orientation  with  the  highest/lowest  score,  the  right/left  2-D  histograms  correspond  to  H+  and  H- . 
(c)  Symmetry  measure  for  a  location  with  a  symmetric  structure  and  a  random  structure. 
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Fig.  4.  (a)-(b)  Dominant  symmetry  direction  and  amplitude  computed  for  multiple  locations  from  a  concertina  wire  and  foliage  data.  In  red  the  most 

likely  symmetric  direction  and  in  blue  the  less  likely,  (c)  Projection  plane  lattice  resolution  influence.  Symmetry  measure  for  one  location,  with  a 
constant  support  region  but  with  different  projection  cell  resolution. 


III.  Approach 

A.  Structure  characterization 

Our  approach  is  inspired  by  the  work  of  Kazhdan  et  al. 
in  [22].  The  authors  propose  to  use  a  reflective  symmetry 
descriptor  to  build  a  signature  of  a  3-D  object  that  can 
be  used  for  evaluating  object  similarity.  Our  approach 
differs  in  a  number  of  areas.  The  authors  deal  with  scenes 
made  of  a  unique  complete  object  (no  missing  data,  no 
clutter).  We  are  dealing  with  partial  data  of  an  object 
of  interest  submerged  in  clutter.  The  authors  deal  with 
object  comparison,  we  are  dealing  with  structure  detection 
and  then  segmentation.  In  addition  we  do  not  compare 
against  a  data  base  model.  Several  methods  have  been 
proposed  to  detect  reflective/bilateral  symmetry  [23],  [24], 
[25],  [26].  We  decided  to  follow  the  approach  of  [8],  [9] 
which  projects  the  data  into  a  different  subspace  to  perform 
comparison. 

B.  Dimension  reduction 

The  core  of  the  method  is  to  project  and  accumulate  the 
3-D  points  of  the  scene  into  a  gridded  planar  surface  to 
produce  a  2-D  histogram  for  each  orientation  of  the  plane. 
Figure  3 -(a)  illustrates  the  geometry  of  the  projection.  For 


each  orientation,  two  2-D  histograms  ( H+  and  //_)  are 
produced  corresponding  to  the  projection  of  the  points  in 
the  two  volumes  ( Vol+  and  Vol-)  separated  by  the  projec¬ 
tion  plane.  The  method  is  characterized  by  the  following 
parameters:  the  width  (W)  and  height  (H)  of  the  planar 
area  where  points  are  projected,  the  maximum  distance  to 
the  plane  for  a  point  to  be  projected  (d^p,  the  number  of 
bins  in  the  projection  plane  (nw  xnh)  and  the  size  of  the 
bins  (sw  xsh),  the  location  of  the  center  of  the  projection 
plane  ( Op ),  and  finally  the  axis  (A)  and  the  angle  (0)  of 
rotation  of  the  plane.  In  this  paper  the  axis  A  is  aligned 
with  the  vertical.  Figure  3-(b)  shows  two  examples  of  2-D 
histograms  computed  for  the  same  location  (center  of  the 
coil)  but  for  different  plane  orientations.  The  projection 
plane  contains  20  x  20  square  bins,  5  cm  wide.  is  set 
to  1  meter.  The  darker  the  pixel  is,  the  higher  the  density 
of  points  in  the  bin  is. 

A  measure  of  bilateral  symmetry  is  naturally  computed 
by  comparing  the  two  2-D  histograms  H+  and  //_.  Several 
distance  measures,  such  as  histogram  intersection  or  the 
Earth  Mover’s  Distance  [27],  can  be  considered.  Here  we 
use  the  normalized  correlation,  denoted  R{H+,H-)  with 
R  G  [—1,1],  which  was  used  by  Johnson  in  [8]  with  similar 


signatures. 

By  computing  this  symmetry  value  for  different  plane 
orientations  we  obtain  a  curve  R  =  /(0).  Figure  3-(c)  shows 
two  examples:  one  from  a  concertina  wire  data  and  the 
other  from  random  data.  The  absolute  maximum  of  the 
curve  determines  the  orientation  and  the  intensity  of  the 
dominant  local  bilateral  symmetry  of  scene. 

In  the  next  section  we  review  the  properties  of  the 
method,  we  present  experimental  results  on  wire  detectabil¬ 
ity,  and  results  on  a  concertina  wire  detection  in  cluttered 
vegetated  terrain. 

C.  Other  approaches  considered 

Other  approaches  were  considered  to  address  this  prob¬ 
lem:  directional  saliency  and  range  image  segmentation. 

1 )  Directional  saliency:  The  idea  is  to  use  not  only  the 
saliency  value  from  the  method  presented  in  Section  II- 
C  but  to  associate  a  direction  with  the  surface-ness  value 
(normal  direction)  and  linear-ness  value  (tangent).  The  two 
critical  issues  are  the  scale  selection  and  the  segmentation. 
Figure  5  presents  two  results  of  the  use  of  directional 
saliency.  In  Figure  Il-C-(a)  a  small  neighborhood  size  is 
selected  and  the  directional  vector  fit  well  to  the  wire 
structure.  In  Figure  Il-C-(b)  a  larger  neighborhood  size 
is  selected  and  the  directional  vector  fit  well  to  the  coil 
structure  rather  than  to  the  wire  structure.  This  example 
is  very  favorable  because  the  ground  has  been  filtered  and 
there  is  no  clutter.  The  second  problem  mentioned  earlier 
is  that  we  are  still  left  with  the  detection  of  the  concertina. 


(b) 


Fig.  5.  Thin  structure  detection  using  directional  saliency.  (a)  Small 
scale  detection  (15  cm  neighborhood),  (b)  Larger  scale  detection  (1.25  m 
neighbohood). 

2)  Range  image:  Obstacle  detection  in  range  images  is 
a  proven  method  for  cross-country  navigation  [3].  Methods 
using  range  image  to  detect  symmetric  pattern  are  currently 
been  considered  because  of  their  potential  for  speeding-up 
the  processing  time. 

IV.  Implementation  and  experimentations 
A.  Properties 

We  now  review  five  properties  of  the  symmetry  measure 
presented  above:  scale,  angular  resolution,  discrimination 
against  noise,  projection  plane  resolution  influence,  and 
finally  robustness  against  clutter. 


1)  Scale:  Symmetry  is  a  scale  dependent  property  of 
objects.  The  scale  is  defined  by  the  size  of  the  projection 
area  (WxH)  and  the  range  of  projection  (d™^).  They  are 
defined  to  capture  the  structure  of  the  object  of  interest.  In 
the  example  presented  in  this  paper  we  choose  d™™-  equal 
to  1  or  1.25  m  and  a  one  meter  square  projection  area. 

2)  Angular  resolution:  In  the  examples  presented  so 
far  the  symmetry  measure  is  computed  for  every  angle  of 
the  rotation  plane,  producing  180  values  per  location.  In 
practice,  the  curves  are  smooth  and  it  is  therefore  possible 
to  evaluate  the  symmetry  for  a  sparse  set  of  orientation, 
thus  speeding  up  the  computation.  In  the  example  presented 
in  Figure  9-(d)  the  angle  interval  is  4°. 

3)  Discrimination  against  random  points:  In  Figure  4, 
we  present  results  showing  that  the  signature  can  discrim¬ 
inate  between  an  organized  points  cloud  from  a  concertina 
wire  and  scattered  points  from  foliage.  For  multiple  lo¬ 
cations,  symmetry  curves  are  computed  and  the  absolute 
extrema  are  extracted.  Figures  4-(a)  and  4-(b)  show  their 
distribution  for  a  concertina  wire  and  foliage  respectively. 
In  Figure  4-(a)  points  are  clustered  around  the  direction 
90°  corresponding  to  the  plane  normal  to  the  main  axis  of 
the  coil.  In  Figure  4-(b)  the  data  is  distributed  randomly. 
Note  the  difference  of  symmetry  intensity  between  the  two 
graphs. 

4)  2-D  histogram  resolution:  Figure  4-(c)  shows  the  in¬ 
fluence  of  the  projection  plane  resolution  on  the  symmetry 
value.  We  use  data  from  a  concertina  wire.  The  range  of 
influence  and  the  size  of  the  projection  area  are  maintained 
constant,  set  both  at  1  m.  Three  pixel  resolutions  are  used 
5,  10  and  25  cm.  As  the  resolution  decreases,  new  modes 
appear  in  the  graph  (118°  and  142°),  the  symmetry  value 
increases  but  the  value  of  the  dominant  symmetry  mode 
remains  at  the  same  location  (90°)  and  its  maximum  well 
above  the  mean  value. 

5)  Clutter:  One  critical  property  for  the  approach  to 
be  applicable  in  vegetated  terrain  is  the  ability  of  the 
signature  to  handle  clutter  around  the  structure  of  interest. 
We  compared  the  symmetry  curve,  for  two  projection 
plane  lattice  resolutions,  for  a  concertina  wire  without,  and 
with  different  levels  of  clutter  simulated  using  a  Poisson 
distribution,  following  [28].  Figure  6-(a)  shows  a  top  view 
of  the  scene  corrupted  with  noise.  The  other  graphs  in 
Figure  6  show  the  symmetry  curve  for  the  same  location, 
for  the  same  range  (1  m),  the  same  projection  plane  size 
(1  m  x  1  m)  but  at  two  different  resolutions  (10x10  and 
20x20).  Note  that,  in  each  graph,  the  curves  are  similar. 
The  scene  contains  2388  points  from  the  concertina.  The 
volume  occupied  is  2.5  m  x  1  m  x  1  m.  The  three  levels 
of  noise  considered  correspond  to  100,  150  and  200  noise 
points  per  cubic  meter. 

B.  Configuration  space  exploration 

In  our  current  implementation,  the  intensity  and  direction 
of  the  bilateral  symmetry  is  computed  on  a  3-D  lattice 
spaced  every  half  to  one  meter.  Figure  7  shows  an  example 
of  such  a  result.  In  this  top  view  of  the  scene  of  an 
isolated  coil  of  wires,  LADAR  points  are  in  red,  the  lattice 


input  data  xyz  Clutter  influence  Clutter  influence 


x  (m)  Angle  (deg)  Angle  (deg) 


(a)  Scene  top  view  (b)  20x20  pixels  2-D  histogram  (c)  10x10  pixels  2-D  histogram 

Fig.  6.  Clutter  influence  on  the  symmetry  measure,  (a)  Scene  with  in  red/blue  the  scene/added  noise  data  (b)-(c)  Symmetry  values  compute  for  the 
same  location  at  low  and  high  resolution.  Each  graph  is  made  of  the  superposition  of  several  plot  corresponding  to  increasing  level  of  noise. 


points  are  in  green  and  the  direction  and  intensity  of  the 
symmetry  is  represented  as  a  blue  vector.  Note  the  change 
of  orientation  of  the  main  axis  of  the  coil  correctly  captured 
by  the  direction  of  the  blue  lines.  The  break  point  is 
located  at  (0.25,3.25).  The  ground  data  points  are  filtered 
out  before  processing. 
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Fig.  7.  Top  view  of  the  scene.  The  red  points  are  the  LADAR  data,  the 
green  points  are  the  location  where  the  symmetry  has  been  evaluated,  in 
blue,  vectors  representing  the  orientation  and  the  intensity  of  the  bilateral 
symmetry. 


C.  Performance  analysis 

1 )  Sensing  issues:  Sensing  thin  structures  is  a  challeng¬ 
ing  problem.  We  tested  several  sensors  against  various  thin 
targets  to  evaluate  wires  detectability.  We  tested  a  Riegl 
LD90  first  pulse  unit,  a  Riegl  LD90  modified  to  extract 
the  full  waveform  [29],  a  LARA  Z+F,  and  a  SICK.  Table 
I  contains  the  characteristics  of  the  wires  tested. 

TABLE  I 

Wires  characteristics 


Color 

Material 

Notes 

Dimention  (mm) 

#1 

Black 

plastic 

flexible 

16  0 

#2 

Copper 

copper 

7  fibers 

5.84  0  (2.032  0) 

#3 

White 

plastic 

9.52  0 

#4 

Grey 

plastic 

flat 

6.98  x  17.90 

#5 

Black 

plastic 

rigid 

12.7  0 

Table  II  contains  for  each  sensor  the  maximum  range 
at  which  two  of  the  wires  are  detectable.  All  but  the  Z+F 
laser  are  Pulsed  laser.  The  Z+F,  a  continuous  wave  laser, 
was  strongly  susceptible  to  mixed  pixels  as  expected. 

TABLE  II 

Maximum  range  detection  in  meter 


Sensor 

Black  wire 

Copper  wire 

Sick 

6.65 

9.00 

Z+F 

* 

* 

Riegl  waveform  sensor 

19.00 

28.00 

Riegl  first  pulse  unit 

23.00 

39.00 

We  used  an  actuated  SICK  LMS-291-S05  mounted  on  a 
cart  to  collect  3-D  scans.  We  performed  a  test  to  determine 
the  detectability  of  a  concertina  wire  using  the  highest 
angular  resolution  setting  of  the  scanner.  With  such  setting 
the  laser  produce  400x800  points  for  a  field  of  view  of 
100°  x  120°.  In  Table  III  is  shows  the  counts  of  range 
measurements  from  the  concertina  wire  scanned  at  different 
distances.  The  ground  was  filtered  and  the  remaining  data 
manually  labeled.  For  some  of  the  same  distances,  Figure 
8  shows  a  front  view  of  the  3-D  data,  with  the  elevation 
color  coded.  Note  that  the  coil  structure  is  preserved  up  to 
8  meters  in  the  configuration  tested.  At  13  m  there  is  only 
spare  returns  from  the  target. 

TABLE  III 

Wire  detectability  by  the  actuated  Sick  laser  at  low 

ANGULAR  RESOLUTION.  THE  RANGE  IS  IN  METER. 


Range 

2.6 

3.6 

4.6 

5.6  6.8 

7.6 

8.6 

9.6 

10.6 

11.6 

12.6 

Points 

2506 

1695 

1178 

820  521 

389 

219 

171 

138 

98 

72 

2)  Scene  segmentation:  We  collected  data  sets  of  the 
coil  positioned  in  a  vegetated  terrain,  seen  from  different 
viewpoints  and  with  different  level  of  clutter.  Figure  9- (a) 
is  an  image  of  the  scene.  Results  from  the  processing  of 
two  data  sets  are  presented  here:  normal  incidence  with 
clutter  (Figure  9-(a/b/c/d))  and  oblique  view  angle  with 
no  clutter  (Figure  9-(e)/(f)).  The  laser  is  positioned  at  7 
m  from  the  coil  and  the  scanner  is  set  in  high  resolution 
mode.  For  space  consideration  we  decided  not  to  present 
the  corresponding  Figure  9-(b/c)  for  the  oblique  view. 
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Fig.  8.  Front  view  of  a  concertina  wire  sensed  at  different  distances 
with  a  Sick  laser. 


as  vegetation.  We  mentioned  this  problem  earlier  in  Section 
II-C.  Figure  9-(d)  shows  the  distribution  of  the  dominant 
symmetry  direction  and  value  for  an  area  surrounding 
the  coil.  The  symmetry  measure  is  calculated  every  four 
degrees.  The  coil  contributes  to  the  accumulation  around 
the  22nd  bin. 

We  produced  similar  results  with  the  oblique  view  of 
the  scene  (Figure-(e/f)).  Figure  9-(e)  is  a  top  view  of  the 
scene.  The  red  points  are  the  LADAR  data,  the  green 
points  are  the  locations  where  the  symmetry  has  been 
evaluated,  in  blue,  vectors  representing  the  orientation  and 
the  intensity  of  the  bilateral  symmetry.  The  distribution  of 
the  same  results  is  presented  9-(f).  The  concertina  wire 
can  be  detected  by  clustering  the  data  in  the  symmetry 
value,  direction  feature  space.  In  this  example,  the  coil  is 
adjacent  to  but  not  cluttered  by  vegetation  unlike  in  the 
previous  example. 

3)  Complexity:  The  current  implementation  is  designed 
for  off-line  processing.  At  this  point  we  consider  all  the 
measurements  collected  from  the  laser,  there  is  no  3- 
D  points  collapse  into  voxels  like  in  [20].  In  addition 
we  sample  the  3-D  space  regularily,  there  is  no  location 
selection  like  during  the  PerceptOR  program.  During  this 
program,  a  similar  approach  of  data  reduction  by  projection 
was  used  to  construct  terrain  signature  for  air-ground 
LADAR  data  registration  for  the  purpose  of  ground  robot 
localization.  The  approach  was  made  more  efficient,  in 
terms  of  computation  for  on-board  robot  application  and 
accuracy,  because  of  a  strategy  for  point  selection.  A  priori 
good  location  canditates  only  were  used  to  compute  the 
signatures  [4]. 


V.  Conclusion 

In  this  paper,  we  investigate  a  method  to  detect  con¬ 
certina  wire  in  LADAR  data.  The  method  allows  the  recov¬ 
ery  of  the  local  dominant  bilateral  symmetry  direction  and 
intensity.  We  show  the  influence  of  the  method  parameters 
using  data  from  a  controlled  environment.  We  also  show 
segmentation  results  of  the  concertina  wire  in  vegetated 
terrain. 

But  the  larger  problem  addressed  here  is  the  one  of 
detecting  organized  thin  structures  in  3-D  LADAR  data. 
This  work  is  a  new  contribution  of  our  effort  started  in  [4] 
and  pursued  in  [20]  to  reliably  segment  vegetation  terrain 
for  ground  mobile  robot  navigation.  Future  work  will 
address  the  characterization  of  additional  gross  structures 
in  addition  to  the  current  bilateral  symmetry.  Statistically 
significant  evaluations  on  various  scenes  are  planned. 


The  data  was  processed  in  three  steps:  1)  ground  sep¬ 
aration,  2)  non-ground  data  classification  in  three  classes 
(clutter,  linear,  surface),  3)  dominant  bilateral  symmetry  de¬ 
tection.  The  ground  data  is  separated  by  using  the  method 
proposed  in  [4]  and  the  result  is  seen  in  Figure  9-(b).  The 
classification  results  of  the  step  2  can  be  seen  in  9-(c). 
Notice  that  it  is  classified  partially  as  linear  and  partially 
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Segmentation  of  non-ground  points  into  linear/scatter/surface  classes  fol¬ 
lowing  [20] 


(b)  Scene  segmented  into  ground/non  ground 


Direction  and  symmetry  value  distribution 


10  20  30  40 

Direction  bin 


(d)  Symmetry  detection 


Segmentation 


-3-2-1  0  1  2  3 

X  (m) 


(e)  Symmetry  detection 


Distribution 


<u 


o> 

E 

E 

CO 


1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 

-0.2 

0  50  100  150  200 


Angle  (deg) 

(f)  Symmetry  value  and  orientation  distribution 


Fig.  9.  Example  of  scene  analysis,  (a)  The  scene:  concertina  wire  against  vegetation  with  branches  added  to  clutter  the  wires,  (b)  Segmentation  of  the 
ground  points  represented  in  light  brown  color.  The  rest  of  the  points  are  colored  based  on  their  elevation,  (c)  Segmentation  of  the  non-ground  points 
into  linear  (blue),  scatter  (green)  and  solid  (red)  structures,  (d)  Symmetry  direction  and  intensity  distribution  computed  at  locations  selected  randomly 
in  the  scene  but  below  one  meter  elevation  from  the  ground.  Normal  incidence  data  collection,  (e)  Top  view  of  the  scene.  The  red  points  are  the 
LADAR  data,  the  green  points  are  the  location  where  the  symmetry  has  been  evaluated,  in  blue,  vectors  representing  the  orientation  and  the  intensity 
of  the  bilateral  symmetry  (f)  Dominant  symmetry  direction  and  amplitude  computed  for  multiple  locations  from  the  concertina  wire  and  foliage  data. 
In  red  the  most  likely  symmetric  direction  and  in  blue  the  less  likely.  (e)/(f)  produced  using  the  side  view  data  collection. 


